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THE LENGTH SPECTRUM OF THE SUB-RIEMANNIAN 

THREE-SPHERE 

DAVID KLAPHECK AND MICHAEL VANVALKENBURGH 


Abstract. We determine the lengths of all closed sub-Riemannian geodesics on the 
three-sphere S 3 . Our methods are elementary and allow us to avoid using explicit for¬ 
mulas for the sub-Riemannian geodesics. 


1. Introduction 

In the case of a compact Riemannian manifold (M, g ) there is a relationship between 
closed geodesics, representing paths of free classical particles in periodic motion, and 
eigenfunctions of the Laplacian A, representing periodic free quantum “waves” (up to 
a phase factor). For this reason, the set of lengths of closed geodesics is called the 
length spectrum , in analogy to the spectrum of the Laplacian. There are in fact precise 
formulas relating lengths to eigenvalues; see for example the announcement of Guillemin 
and Weinstein for a readable discussion with references [12], 


So far there is no such formula relating lengths and eigenvalues in the case of a com¬ 
pact sub-Riemannian (sR) manifold. We recall that a sR manifold is a manifold with a 
specified linear subbundle hi (the “horizontal bundle”) of its tangent bundle, along with a 
Riemannian metric on hi. Distances between points are then measured using curves that 
are constrained to have tangent vectors in hi (“horizontal curves”). In fact, when hi is 
the span of a set of bracket-generating vector fields, then the Chow-Rashevskii theorem 
says that any two points are connected by a horizontal curve, a result that even experts 
find surprising (m, p.178); thus given any two points there is a shortest horizontal curve 
connecting them; it is called a sR geodesic. 


sR geometry is of practical interest; for example, the problem of parallel parking a car, 
or, even worse, a car with a trailer, is a problem in sR geometry [Tj.|16j. And there are 
further surprises from the purely mathematical point of view, one being Montgomery’s 
proof of existence of singular sR geodesics (singular in the sense that they do not satisfy 
the geodesic equations (Hamilton’s equations)) |1]1],|T5]. This and other relatively recent 
results in sR geometry then inspire renewed interest in the subLaplacian: the operator 
naturally associated with the given (sub-)Riemannian metric on hi. 


In this paper, with the goal of understanding a single example, we compute the sR 
length spectrum of the three-dimensional sphere S 3 with its standard sR structure; this is 
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to be compared with the spectrum of the subLaplacian on S 3 , known by Taylor |19] and 
generalized to other connected, semisimple Lie groups by Domokos [D]. We expect that 
a general theory relating the sR length spectrum to the spectrum of the subLaplacian 
would be amenable to the tools of microlocal analysis as in the Riemannian setting; the 
work of Cohn de Verdiere, et ah, gives hope that this will be accomplished [6]. 


We focus on S 3 with its standard sR structure because it is perhaps the simplest 
compact manifold with a sR structure, and there are no singular sR geodesics on S' 3 ; that 
is, all sR geodesics arise as projections of solutions of Hamilton’s equations [T5j. Moreover, 
we wish to compare the sR setting to the Riemannian setting, in which the spheres S n 
are of fundamental importance, as examples of manifolds all of whose geodesics are closed 
and have the same length T ; in general this is equivalent to most of the spectrum of y —A 
being concentrated near an arithmetic progression =f-k + f3, k — 1, 2,..., for some constant 
/3 [TD]. As we will see, in the case of S 3 not all sR geodesics are closed, and not all have 
the same length: 

Theorem. The set of lengths of the closed sR geodesics on S 3 is 

{27 Ty/n; n G N}. 


Others have studied the sR geodesics on S 3 [2J, [Tj, [13] (see also the survey article 0), 
but we compute their lengths and differ from the previous work in that we consistently 
use Hopf coordinates on S 3 and avoid using explicit formulas for the sR geodesics; we 
believe it clarifies the presentation to not use explicit formulas. 


We introduce the sR structure and geodesic equations in Section [2] using Hopf coordi¬ 
nates, and in Section [3] we categorize the qualitatively different types of sR geodesics. In 
Section [I] we determine which sR geodesics are closed, and in Section Owe compute their 
lengths, resulting in the theorem above. Finally, in Section [6] we compare the sR length 
spectrum to the previously-known spectrum of the subLaplacian. 


This work was supported by a SURE (Summer Undergraduate Research Experience) 
Award at California State University, Sacramento. 


Remark. During peer review, it was pointed out that the above result is contained in 
the paper [5] (see their Theorem 2). However, our proof is entirely new and has the 
advantage of being elementary after the introduction of Hamilton’s equations (]2]) in our 
chosen coordinate system. 
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2. S 3 in Euclidean and Hopf Coordinates 
First we consider S 3 as a subset of M 4 : 

S 3 = (Oi, Vi, X2, 2/2) e K 4 ; x\ + y\ + x\ + y\ = 1 }. 

On S 3 we have the orthonormal vector fields 

0 0 0 0 

» : = “2/iu--2/2 o- 

<7£i <7?/i (7X2 <7£/ 2 

^ 0 0 0 0 
E\ := —X2~^ - V 2 / 27 ;-6 ^it;-2/iti— 

<7Xi (7?/i OX 2 01/2 

0 0 0 0 

E 2 ■= 2/2 Tj-2127^-h 2/l~ -1- IlJ- 

< 7 x 2 ui)2 

which satisfy the Lie bracket relations 

[V, #i] = —2£ 2 , [£ 2 , H] = -2 E u [Ei, E 2 ] = -2V. 

Thus T-L (S' 3 ) = Spanjiy, E 2 ) is a bracket-generating tangent subbundle, and by the Chow- 
Rashevskii Theorem any two points on S 3 are connected by a sR geodesic. 


The orbits of the flow generated by V are the circles of the Hopf fibration [3], so we 
find it convenient to use Hopf coordinated on S 3 : 

Xi = cos 9i sin 9 0 
2 /i = sin 9i sin 9 0 
x 2 = cos 9 2 cos 6*o 
y 2 = sin 9 2 cos9 0 , 

for 0 < 9 0 < t/ 2 and 0 < 9j < 2 tt , j = 1,2. We picture the (6b, Qi, 6 l 2 )-space as “the 
Hopf cube” (0, 7 t/2) x (0, 27t) x (0, 2i r). When we have occasion to exit the Hopf cube, we 
simply return to the definition of Hopf coordinates to make the correct interpretation: 

(i) For the 9\ and 9 2 coordinates the values 0 and 27 t are identified. 

(ii) When a point crosses the 0o = 0 plane we have that 6b changes direction (“bounces”) 
and (6*i, 0 2 ) is identified with ( 9 X + 7r, 9 2 ). 

(iii) When a point crosses the 6b = 7 t/ 2 plane we have that 9q changes direction and 
(9i, 9 2 ) is identified with (0 1; 9 2 + tt). 


The (round) Riemannian metric in Hopf coordinates is: 

(1) ds 2 = d9l + sin 2 6 l 0 d9\ + cos 2 6b d# 2 , 


and the Laplacian is 


A 


1 0 
sin (2 9 0 ) 09o 


o 


0 d 2 

sin(26» 0 )— + csc 2 6» 0 + sec 2 6> 

(7C7o uU]_ 


9 2 

w 


* “3-sphere.” Wikipedia, The Free Encyclopedia. Wikimedia Foundation, Inc. June 16, 2015. Web. 
July 1, 2015. 
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We now write the sR structure in Hopf coorcli 
coordinates to M 4 , allowing us to write the -JE 

to functions on S 3 we get: 

Odd 

—— = cos 6\ cos 0 o — - sin 6\ esc 0 0 — 

OX i OUq OU] 


in terms of the 

oQj 


'O 

9 9 9 

—— = sin cos 0 o —— + cos 0i esc 0o 
<9r/i Wo 


90 


'i 


9 


Our vector fields are then 

v = — + — 

90! 902 


■S— = - cos 0 2 sin 00 7^- - sin 0 2 sec 0 O7 ^- 
(7X2 (7C7o Uv2 

d d d 

—— = - sin 02 sin 0 o — + cos 0 2 sec 0 ot^- 
oy 2 Ouq OU 2 


Ex 

E- 


UUi UU2 

d d 

= - cos(0i + 02 ) 7 ^ 7 - + sin(0! + 0 2 ) cot 0 O —-sin(0! + 0 2 ) tan0 o 

(7C7o C/t7l 

9 


_d_ 

'0 oth 90 2 

d c) c) 

sm(0i + 02 )ttt- - cos(0i + 0 2 ) cot 0 O — + cos(0i + 0 2 ) tan0 o —. 

C/uo Ou\ OV2 


E 2 — ~ 8111(0! + 02) 7:7 -COS(01 + 02) cot 0 O — + COS(0i + 0 2 J tai 

90o 90i 

The commutation relations hold, same as before, and the vector fields are still orthonormal 
(of course, with respect to the Riemannian metric in Hopf coordinates). 


The sR metric, written in Hopf coordinates, is 

(10 0 \ 

S — I 0 cos 2 0 o sin 2 0 o — cos 2 0 o sin 2 0 o J . 
yO — cos 2 0 o sin 2 0 o cos 2 0 o sin 2 0 o J 

Indeed it is easy to check that E\ and E 2 are orthonormal with respect to S, and V is in 
the kernel of S. Written as a two-tensor, 

S = 00 o ® 00 0 + cos 2 0 o sin 2 0 o (00 1 — 00 2 ) ® (00i — 00 2 ). 


The sR Laplacian, written in Hopf coordinates, is 


A-sR — E\ + E 2 — 


sin(20 o 


.J_ osin(2()o )J_ + 


9 9 

cotflo^-tanOo — 


We can consider the sR metric as being the limit of certain penalty metrics, where the 
V direction is penalized by a factor A > 1. After simple linear algebra (multiplying the V 
direction by A, multiplying the other directions by 1, and then applying the Riemannian 
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is given 


metric), the A-penalty metric 

/l 0 


Px = 


0 

cos 2 #!) sin 2 0 o 


\ 


0 (A 2 - 1) sin 4 + sin 2 0 o (A 2 - 1) 

^0 (A 2 — 1) cos 2 0 o sin 2 0 o (A 2 — 1) cos 4 0 o + cos 2 0 o ) 

Indeed, one can check that in fact V, E \, and E 2 are orthogonal with respect to this 
metric, that E 1 and E 2 have length 1, and that V has length A. We can easily compute 

det P\ = A 2 cos 2 0 o sin 2 0 o 


and 


/I 


P^ = 


0 


0 \ 

A -2 — 1 


0 cot 2 0 o + A 2 
^0 A” 2 — 1 tan 2 0 o + A _2 y 


From this we find that the A-penalty Laplacian on S 3 is: 

M2eo) w 0 + ( coce °ik ~ + A ~ 2 (Jk + w 2 . 


d 2 


2 cot 


Aa d0 2 

That is, 

A a = E 2 + El + A" 2 !/ 2 , 

might have been expected. 


as 


Montgomery discovered an example in which geodesics with respect to the A-penalty 
metric converge (as A —> oo) to sR geodesics that do not solve the sR geodesic equations, 
in contrast to the Riemannian setting; that is, Montgomery discovered so-called singular 
geodesics [23 . For the case of S 3 (and more generally, in the contact case), singular 
geodesics do not exist, so it suffices to study the geodesic equations, or, equivalently, 
Hamilton’s equations [IB] . 


We denote the dual variable to 9j by The sR Hamiltonian is then 

H(6,0 = |£o + |(cot0 o 6 - tan0 o £ 2 ) 2 - 
Hamilton’s equations, giving the sR geodesics, are then, for j = 0,1,2, 

dH 


9j 


Zj = 


3 

dH 


DO* 
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do = £o 

d\ = cot 2 6*0 fi - 6 
e 2 = tan 2 6*0 £ 2 - 6 
4 = cot6 0 csc 2 6* 0 £ 2 — tan0 o sec 2 6* 0 £ 2 

6 = o 
6 = o. 

One obvious advantage of using Hopf coordinates is that 6 and £ 2 are constant along 
the flow; in addition, as always H is constant along the flow, so we already have three 
conserved quantities. Also, these equations have a clear symmetry; for example, 

cot(| — 6q) csc 2 (| — 6 0 ) = tan 6 *o sec 2 6 0 . 

The penalty Hamiltonian is: 

H X ( 9 ,£) = H + 5^(6+ 6) 2 

= + |(cot 6* 0 6 — tan # 0 6) 2 + 2^2 (6 + 6 ) 2 - 

The corresponding penalty Hamiltonian equations, giving the penalty geodesics, are then 

#0 = 6 

0 i = cot 2 6 *o 6-6 + A- 2 (a + 6 ) 

62 — tan 2 6*0 £2 — £1 + A 2 (6 + 6) 

(4) 

£0 = cot 6* 0 esc 2 d 0 £ 2 — tan 6*0 sec 2 6* 0 £ 2 

6 = 0 
6 = 0. 

For the case of the Riemannian metric on S 3 , that is, the case A = 1, the equations 
simplify, and we get 

0\ = esc 2 0 o £1 

0 2 = sec 2 0 o 6 - 

When A = 1, the solutions of Hamilton’s equations are great circles on S 3 . 


Explicitly, 


( 2 ) 


3. Categorizing sR Geodesics 

Our categorization of sR geodesics is based on a reduced problem. In Hamilton’s 
equations ([2D, since 6 and £ 2 are constant along the flow, we can isolate the equations 

$0 = 6 

£0 = cot 6* 0 csc 2 6* 0 £ 2 — tan 6*0 sec 2 d 0 £2 > 
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which are Hamilton’s equations for the sR Hamiltonian H considered as a function of two 
variables 

(5) H{6 o,f 0 ) = + |(cot6»o6 - tan O 0 t, 2 ) 2 - 

Equation (J5J) can be viewed as a one-dimensional energy equation: it is of the form 

Energy = Kinetic Energy + Potential Energy, 
with potential function 

U = |(cot# 0 £i — tan^o^) 2 - 
We now list the various disjoint cases: 


1. A fixed point in the (#o,£o) phase plane. From our original choice of coordinates 
we may assume that 6 0 = 7 t/4 , and then 

la. £i = £ 2 . (This is precisely the case when H = 0.) Then from Hamilton’s 
equations 9 0 , 9\ , and 0 2 are constant; this gives a degenerate sR geodesic of 
length 0. 

lb. = —£ 2 0. Hamilton’s equations then say that the speed on the Hopf 

cube is \/2 |£i — £ 2 |, and the length of the (simple) closed curve on the Hopf 
cube is a /2 27 t, so the period is 27 t /|£ 1 — £ 2 |. On S 3 the speed is — £ 2 |, so 
the length of this closed sR geodesic is 27 t. See Figure [D 


We categorize the remaining cases in terms of the potential function U. 


2. The “free” case U = 0. This happens precisely when £i = £2 = 0. (We have 
already dispensed with the case when 9 0 is constant.) By Hamilton’s equations 
9i, 0 2 , and £ 0 ) are also identically zero, while 9 0 = £ 0 - That is, we have a point 
with speed |£ 0 | moving purely in the ^-direction; the length of this (simple) closed 
geodesic is 27 t. (It is both a geodesic and a sR geodesic.) See Figure [2j 

3- £1 7^ 0 and £ 2 7^ 0. Then U is a potential well with a single non-degenerate mini¬ 
mum occurring when tan 4 d 0 = Typical potential functions are in Figures [3] 

and[4j for £1 and £ 2 with the same and opposite signs, respectively. 


Since in this case 0q is not constant, its period is 

rb d0 n 


Period(6 l o) = 2 


= 2 


. V 2 ( H - u ) 

* dd „ 


L \J2H - (cot9 0 £1 - tan0 o £ 2 ) 2 
Here a and b are the “turning points,” where the kinetic energy is zero. 
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Figure 1. Case lb. Figure 2. Case 2. 




Fortunately it is possible to evaluate this integral using freshman calculus. Sub¬ 
stituting 


x = cos 2 6* 0 , 0 < 6 0 < 7r/2, 


we get 


Period(# 0 ) 


r»cos 2 (a) 


dx 


cos 2 (6) 


— (£i + £ 2 ) 2 ]£ 2 + 2 (H + + £,%)x — £2 


The limits of integration are exactly the points where the denominator vanishes 
(where the velocity is zero), and we recall that the Hamiltonian for Riemannian 
geodesics is H\ — H + + £ 2 ) 2 (the case A = 1), so we have 


i r os2(a) dx 

J cos 2 (6 ) ^/(cos 2 (a) — x)(x — cos 2 (b)) 


Period(6 , 0 ) 
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Figure 5. Case 3. Ex. 1. 


Figure 6. Case 3. Ex. 2. 


This is an integral known to be solvable by elementary functions. Following page 
366 of Advanced Calculus (New Edition) by Frederick S. Wood^ [21], we make the 
substitution defined by 


cos 2 (a) — cos 2 ( 6 ) 


^ + 1 — -UUT 

x — cos z (b) 


and finally get the answer 



7T 




For future reference, we note that this is one-half the period of the (Riemannian) 
geodesic flow; after all, the speed of the geodesic flow is y/2H[, and we know the 
length of each geodesic, a great circle on S 3 , to be 27 t. (See Section HI) 

Examples are pictured in Figures [5] and [ 6 ] where £| and £ 2 have the same and 
opposite signs, respectively. 

4. It remains to check the exceptional cases when (£i = 0 and (2 7 ^ 0} and when 
{£1 7 ^ 0 and £2 = 0}. For example, when £1 = 0 the potential function is 

U = | tan 2 6*0 £f, 0 < 9 0 < tt/2. 

fThis is the book Feynman mentions in Surely You ’re Joking, Mr. Feynman! as giving him valuable 
tricks for integration m- 
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Figure 7. 0 < £i <c £ 2 . Figure 8. 0 < £ 2 < £i- 


The force induced by this potential causes the point to exit the Hopf cube through 
the do = 0 plane; rather we interpret it as bouncing off the plane, returning to the 
Hopf cube but with 6\ shifted by tt. (See Section El) With reasoning as in the 
previous case, we find that again Period(#o) = ir/y/2Hi. The case when (i ^ 0 
and £2 = 0 follows by renaming the variables 9q -h- f — &o and £1 ^ £ 2 . In Figures[7| 
and [ 8 ] we have the cases when 0 < £1 £2 and 0 < £2 -C £1, respectively, which 

illustrate how exiting the Hopf cube and re-entering after a 7 r-shift appears as a 
limiting case. 

Finally, we note that all sR geodesics are simple curves; that is, they do not self-intersect 
except trivially for closed curves. In Cases 1 and 2 above it is obvious. In Cases 3 and 4 
we only need to wait until £0 is zero, corresponding to the d 0 -particle having zero kinetic 
energy in the potential well U. When (60,61,62) returns to that value, clearly £ 0 is zero 
again, £i ,£ 2 are the same as always, and the 9j and £,■ return to their values; thus the 
curve only self-intersects in the case of a closed curve, at the end of a period. 


4. Determining Which sR Geodesics are Closed 

In this section we identify the closed sR geodesics on 5 3 ; we only need to consider the 
Cases 3 and 4, and we may assume that the initial value of £ 0 is zero. (See the comment at 
the end of Section El) Hurtado and Rosales P3J found a necessary and sufficient condition 
in terms of geodesic curvature (see also the expository paper of D’Angelo and Tyson my- 

Theorem. |13] Let 7 : M —> S ' 3 be a complete sR geodesic of curvature X. Then 7 is 
a closed curve diffeomorphic to a circle if and only if A/\/l + X 2 is a rational number. 
Otherwise 7 is diffeomorphic to M and is dense in some group translate of a Clifford torus. 
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Their proof relies on closed-form expressions of the sR geodesics. Here we give a 
condition which does not rely on closed-form expressions. 


From the A-penalty Hamilton’s equations (J3J), we see that the sR Hamiltonian vector 
field for the Hamiltonian H is the difference of the Hamiltonian vector fields for the 
Hamiltonians Hi and Hy = |(6 + 6) 2 - Moreover, the vector fields Lie-commute (it is 
easy to see that the Poisson bracket of H\ and Hy is zero), so the Hamiltonian flows for 
H\ and Hy commute. We can thus consider the H- flow as an LR-flow followed by an 
Hy-RoW. 


The Hamiltonian for the Riemannian geodesics may be written as 

Hi(6, f) = + |(csc 2 0 o + sec 2 6 0 £ 2 ) , 

(the penalty Hamiltonian (J3J) with A = 1), so the first of Hamilton’s equations, giving the 
velocities, are then 

8o = 6 
9i = csc 2 8 0 6 
8 2 = sec 2 0 o 6- 

We see that the speed (measured using the Riemannian metric (JT]) ) is \/2II\. which is 
constant. Moreover, the length of the Riemannian geodesic is 27T, being a great circle, so 
that the period of the closed orbit is 2ir/y/2Hi = 2 x Period(6 , o). 


On the other hand, the Hamiltonian Hy has Hamiltonian equations 

6 = 6 + 6 
02 = 6+6 
6 = 6 = 0. 

Thus the speed (with respect to the Euclidean metric on the Hopf cube) is \/2|6 + 61 - 
The length of the orbit (a circle fiber of the Hopf hbration) is \/2- 2n r, so the period of the 
Hy- flow is 27 t/| 6 + 61- [If might seem strange that we find the speed and length with 
respect to the Euclidean metric on the Hopf cube, but the Euclidean metric is sufficient 
to compute the period of the Hy- flow.] 


For a combination of an H \-flow and an Hy-Row to result in a closed curve, we need the 
HyRow to return Oq to its original value (since the Hy- flow has no component). Thus 

the time elapsed must be an integer multiple of Period(#o) = 'Kj\/2Hi. If the integer is 
odd, the Hy -flow takes the point to its antipodal point, and we would need a half-period 
of the Hy-Row to return to the starting point. If the integer is even, the Hy flow takes the 
point back to itself, and we could only allow full periods of the Hy- flow. To summarize, 
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a necessary and sufficient condition for a closed sR geodesic is: 


11 7T 7T 

time elapsed = p x — -— = q x 


16 + 61 

where p, q G {1, 2, 3,...} are either both odd or both even. In particular, 

P 16 + 61 


( 6 ) 


__ = \le Qn ( 0 , 1 ), 

q y/2H[ V H x V ^ 


The quantity p/q is conserved along the flow and is positively homogeneous of degree zero 
in the 6 var iables. The condition © is also sufficient to have a closed sR geodesic. If it 
holds, then we have 

7T 7r 

H -Period = p x --- = q x 


16 + 61 11 " y/m' 

for the least such integers 0 < p < q that are either both odd or both even. 


When plotting sR geodesics in Cases 3 and 4, we can fix any r e Qfl (0,1) and re-write 
the closure condition (JBJ) as 


f 2 - 
so — 


(6 + 6) s 


— esc 2 d 0 — sec 2 6*0 


We can always hnd initial conditions satisfying this. Indeed, in Case 3 we can take any 
nonzero 6 and 6 and then take 6q to maximize the right-hand side: tan 2 $o = 16/61- If 
6 and 6 have the same sign, the right-hand side is always positive. If 6 an d 6 have 
opposite signs, we need 


6 + 6 

6-6 


> r, 


which is only valid for certain 6 and 6 - Case 4 is similar. Then we can solve for £ 0) nse 
those numbers as the initial conditions in Hamilton’s equations, and then plot the closed 
sR geodesic. Taking, for example, r = 1/5, 6 — 0.6, and 6 — 0.7 we get the sR geodesic 
in Figure EB 


5. The sR Length Spectrum 

To calculate the lengths of the closed sR geodesics we again only need to consider 
Cases 3 and 4 (the cases where 9 0 oscillates). We found in the previous section that a 
sR geodesic is closed when the period of the Hi -flow and the period of the Hy-Q ow are 
commensurable. Then we have 


7T 7T 

16 + 61 = qX Vzh~i 


(7) 


Period of H -flow = p x 
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Figure 9. An example with r — 1/5. 


for the least such integers 0 < p < q where p, q are either both odd or both even. Since 
we know the speed of the sR geodesic is a constant V“2H, we have that the length is 

Length = Period x Speed 


irq 


( 8 ) 


x/2 Hi 


x 


V 2 H 


H 





p 2 , 


for the least integers 0 < p < q satisfying (J7J) where p, q are either both odd or both even. 


We have another formulation of length that explains the repeating patterns seen in the 
figures. We know that the distance traveled in one do-period is 

Period(d 0 ) x speed = tt 

Thus the length of a closed sR geodesic is 

Length = 7r x (number of do-oscillations) x 

Comparing with equation (jHJ), we find that 

number of do-oscillations = q. 

Moreover, we see from Hamilton’s equations that the curve segments traced out by do- 
oscillations are congruent to each other. A similar argument shows that Riemannian 
geodesics in Hopf coordinates consist of two d 0 -oscillations, as illustrated in Figure flOl 
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Figure 10. A Riemannian geodesic in Hopf coordinates. 


To summarize, we have found that if a sR geodesic is closed then the initial conditions 
must satisfy 


and that the length of the closed sR geodesic is 

Length = 7 t\J q 2 — p 2 , 

for the least integers 0 < p < q satisfying (JTJ) where p, q are either both odd or both even. 


In fact, every such number is attained as a length; we simply follow the procedure: 

(i) Choose any p/q G Q fl (0,1), with gcd(p, q) = 1. 

(ii) As seen at the end of Section 01 we can choose initial conditions so that 

p = r w i6+6i 

q V H\ ^2H[ ' 

Thus 

7T 7r 

V X 77-~T = q X ; 

16 + 61 V2Th 


(iii) If p and q are both odd, then the sR geodesic with those initial conditions has length 
7T \Jq 2 — p 2 . If one of {p, q} is odd and the other is even, the sR geodesic with those 
initial conditions has length 2ir^/q 2 — p 2 . 
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Thus the length spectrum consists of 2n and the numbers 

TT\/q 2 - p 2 

where 0 < p < q are odd integers with gcd(p, q) — 1, and 

2ta J q 2 — p 2 

where 0 < p < q are integers, one odd and the other even, with gcd(p, q) — 1. 


We now give an alternative characterization of these numbers. It is simpler to work 
with squares of lengths divided by 7r 2 . Then we wish to characterize the set S of numbers 
consisting of 4 and 



where 0 < p < q are integers with gcd(p, q) — 1 and 


e 


1 if p and q are both odd 
4 if one of {p, q] is odd and the other is even. 


In the e = 1 case we take the examples p = 2k — 1 and q — 2k + 1, k G N, to get 

q 2 — p 2 = 4(2fc), k G N. 

In the e = 4 case, we take the examples p = k and q — k + 1, k G N, to get 

4(g 2 -p 2 ) =4(2fc + l), keN. 

This shows that 4N C S. Now suppose that n e S and 4 { n. Then clearly n can only be 
in the e = 1 case, so there would be odd integers 0 < p < q with gcd(p, q) — 1 such that 
n = q 2 — p 2 . This is easily seen to be impossible. Thus in fact 4N = S. 


We note that if n G S and 8 | n, then n cannot be in the e = 4 case, and that if n e S 
and n = 4(2 k + 1), k G N, then n cannot be in the e = 1 case. Both of these statements 
easily follow from parity arguments. 


Converting back to the language of lengths, we find that the set of lengths of the closed 
sR geodesics is 

{27r\/n; n G N}. 

By the previous paragraph, odd n correspond to “full periods” of the H v -flow and geodesic 
flow (the e = 4 case), and even n correspond to “half periods” of both the H v -flow and 
geodesic flow (the e = 1 case). 
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6. The Spectrum of the subLaplacian 


The subLaplacian — A sR has a compact resolvent, hence has a pure discrete spectrum 
0 < Ai < A 2 < • • • < \ n < ■ ■ ■, with A n — > Too as n —* Too, and a complete orthonormal 
set of eigenfunctions. (See, for example, the recent paper of Colin de Verdiere, et ah [hj.) 
In fact, in the case of S 3 , the eigenfunctions of the subLaplacian are the same as the 
eigenfunctions of the Laplacian. We recall that Aa = Ef + E% + A~ 2 V 2 is the A-penalty 
Laplacian, with A = 1 giving the Riemannian Laplacian on the sphere A S 3, and A = oo 
giving the subLaplacian on the sphere. In Hopf coordinates we have 


A sR — E 2 + E\ — 


sin( 2 d 0 ) d 6 ( 


■ o smi/tAd 




It is easy to see that V = commutes with A sR , hence A sR commutes with A S 3 . 

Thus A sR and A 53 have a common complete orthonormal set of eigenfunctions am 
the eigenfunctions of A sR are simply the spherical harmonics. 


Particularly noteworthy is (x\ + iyi) k = s'm k 9 0 e lkei . It is a “Gaussian beam”: a family 
of eigenfunctions of both A 33 and A sR that concentrates along a great circle. Zelditch 
singles out this example in the Riemannian setting (see, e.g., p.76). It would be 

interesting to see if it is possible to construct, localized to each sR geodesic, a quasimode 
or Gaussian beam in the spirit of Ralston nu, pi. 


Taylor used the Peter-Weyl Theorem to find the eigenvalues of A sR (T9]; Domokos gen¬ 
eralized, using subelliptic Peter-Weyl and Plancherel Theorems on compact, connected, 
semisimple Lie groups [9]. To summarize, the eigenvalues of — A 53 are m{m + 2) for 
m G {0,1, 2,...}, and the eigenvalues of —A sR are (for the same m; the operators have 
the same complete orthonormal set of eigenfunctions): 

4 mj — Aj 2 + 2m, j G {0,1, 2,..., m}. 

For reference, the eigenvalues of the A-penalty Laplacian 

—Aa = —A sR — A 2 V 2 
are 

(1 — A~ 2 )4j(m — j ) + m(2 + A ~ 2 m), 
for m G {0,1,2,...} and j G {0,1, 2,..., m}. 


At this point we will not conjecture a general formula relating the sR length spectrum 
of a bracket-generating compact sR manifold (which for S 3 is {2i Ti/r; n G N}) to the set 
of eigenvalues of the subLaplacian counted with or without multiplicities (which for S 3 is 
{ 2 m; m = 0 , 1 , 2 ,...}). 
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